Measurement of the 8 B Solar Neutrino Flux with the KamLAND Liquid Scintillator Detector 
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We report a measurement of the neutrino-electron elastic scattering rate from 8 B solar neutrinos based 
on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 
5.5 MeV analysis threshold is 1.49±0.14(stat)±0.17(syst) events per kton-day. Interpreted as due to a 
pure electron flavor flux with a 8 B neutrino spectrum, this corresponds to a spectrum integrated flux of 
2.77±0.26(stat)±0.32(syst) xl0 6 cm _2 s _1 . The analysis threshold is driven by 208 T1 present in the liquid 
scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 1 Be. The 
measured rate is consistent with existing measurements and with Standard Solar Model predictions which in- 
clude matter enhanced neutrino oscillation. 

PACS numbers: 26.65.+t,14.60.Pq,13.15.+g 



Neutrinos from 8 B beta-decay (Q~18MeV) dominate the 
high-energy portion of the solar neutrino spectrum and have 
been a crucial source for solar neutrino experiments. Such 
experiments, together with reactor anti-neutrino experiments, 
have established flavor change through Large Mixing An- 
gle (LMA) neutrino-oscillation with matter effects, introduced 
by Mikheev, Smirnov and Wolfenstein (MSW) JH 0], as a 
consistent resolution of the solar neutrino problem 1131—1 1 1 - 
Measurements of neutrino-electron scattering using water- 
Cherenkov techniques have been achieved in a number of dif- 
ferent detectors: Kamiokande [ U , 12 ] , Super-Kamiokande JH 
H0] and SNO JHi3l- An advantage of this detection tech- 
nique is sensitivity to the neutrino direction of incidence, al- 
though analysis thresholds of ~5 MeV are typically imposed. 
Recently SNO has achieved a threshold of ~3.5MeV lioll . 
Liquid scintillator detectors may also be used to measure 



neutrino-electron elastic scattering; while insensitive to the 
neutrino direction, these detectors have better energy resolu- 
tion and have the possibility of lower energy thresholds. A 
measurement with a threshold of ~ 3.0 MeV, the lowest to 
date, has been reported by Borexino IU4II . In this paper we 
report an independent, liquid scintillator based measurement 
of the 8 B flux with KamLAND. 

The KamLAND detector consists of 1 kton of liquid scintil- 
lator (LS) confined in a 6.5-m-radius, 135/im-thick, spherical 
balloon. The LS is 80% dodecane, 20% pseudocumene, and 
1.36±0.03g/l of PPO; the density is 0.78 g/cm 3 . The balloon 
is suspended in purified mineral oil within a 9-m-radius stain- 
less steel sphere (SSS). The scintillation light is recorded by 
an array of 1879 photomultiplier tubes (PMTs) mounted on 
the inner surface of the SSS; 554 are reused 20-inch PMTs 
from Kamiokande, the remainder are new 20-inch PMTs with 
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TABLE I: Production rate of muon-spallation isotopes that dominate 
the background (in events per kton-day) lfl5ll . Non-bright muons 
produce less than 7 x 10 5 p.e. in the PMT array. 



Isotope 


Q[MeV] 


">"l/2 [S] 


All Muons 


Non-Bright Muons 


8 Li 


16 


0.84 


15.6±3.2 


0.7±0.4 


8 B 


18 


0.77 


10.7±2.9 


<0.006 


"Be 


11.5 


13.8 


1.4±0.3 


0.27±0.30 



the photo-cathode masked down to 17-inch. The 17-inch 
tubes have better timing and single-photo-electron (spe) reso- 
lution. A 3.2 kton, water-Cherenkov detector surrounds the 
SSS. It is used as a muon anti-coincidence counter and as 
shielding against external 7-rays and neutrons. 

Event position (vertex) and energy are reconstructed based 
on the photon hit- time and charge distributions and a detec- 
tor response model which is calibrated by periodically de- 
ploying radioactive sources: 203 Hg, 137 Cs, 68 Ge, 65 Zn, 60 Co, 
241 Am 9 Be, 210 Po 13 C. For the current analysis, the recon- 
struction omits the 20-inch PMTs due to their poorer tim- 
ing and spe resolution; in this mode the energy resolution is 
7%/ v /E(MeV). 

The signal of interest, neutrino-electron elastic scattering 
from 8 B solar neutrinos, produces recoil electrons with recon- 
structed energy up to ~20 Me V. We construct a vertex-^ 2 test 
to select electron-recoil-like events. The test compares the ob- 
served PMT charge and hit time distributions as a function of 
the distance to the event to the expected distributions. The ex- 
pected distributions and selection criteria were developed and 
calibrated based on source calibration data and muon-tagged 
/3~events from cosmogenic 12 B/ 12 N . The efficiency of the 
selection is 0.99 ± 0.01. 

The analysis is carried out in reconstructed energy with a 
threshold of 5.5 MeV, which corresponds to a 5 MeV thresh- 
old in physical electron recoil energy. The threshold is chosen 
to reject background from 208 Tl(Q=5.0MeV, T!/ 2 =3.05 min) 
produced by residual 232 Th contamination present in the LS 
at a level of (7.90±0.25) x 10~ 17 g/g. 

An energy-scale model was developed to convert physical 
energy to reconstructed energy. The model includes the lin- 
ear response of the LS and particle-dependent non-linearities 
such as scintillator quenching and Cherenkov light produc- 
tion. The energy scale is constrained by calibration data, 
tagged a-particles from 214 Bi- 214 Po/ 212 Bi- 212 Po decays from 
residual 238 u/ 232 Th in the LS and high energy f3~ decays (Q 
~ 14 MeV) from cosmogenic 12 B/ 12 N . 

The measurement reported here is restricted to data ac- 
quired before the KamLAND LS purification. The data was 
taken between April 2002 through April 2007, and corre- 
sponds to 1432.1 days of run-time. To reduce the external 7- 
ray background the fiducial volume (FV) is reduced to a 6-m- 
high, 3-m-radius cylinder centered at the detector center. This 
shape defines an optimal analysis volume that takes advantage 
of further shielding by the ultra-pure LS from backgrounds 
coming from the cavern's rock walls and a large steel deck 



TABLE II: Estimated background contributions for the full exposure 
and after all cuts. The shorter lived spallation products 12 B, 12 N, 9 C, 
9 Li, and 8 He are also considered and their contributions are summa- 
rized by "Spallation Other". 



Background 


Counts 


Spallation n Be 


89.1 


+ 


19.1 


Spallation 8 Li 


20.5 


± 


4.0 


Spallation 8 B 


11.0 


± 


3.0 


Spallation Other 


0.4 


+ 


0.6 


External gamma rays 


25.2 


± 


12.6 


8 B CC on 13 C GND 


5.8 


± 


1.4 


Reactor 77 e 


1.6 


± 


0.1 


8 B CC on 13 C 3.51 MeV 


1.1 


± 


0.4 


hep ES 


0.6 


± 


0.1 


Atmospheric v 


2.0 


± 


2.0 


Total 


157.3 


± 


23.6 



which caps the experiment. The FV fraction and systematic 
uncertainty is determined from short-lived muon spallation 
products, mainly 12 B(Q=13.4MeV, t 1/2 =20.2ms), which 
are assumed to be produced uniformly in the LS. The FV frac- 
tion is taken as the ratio of the number of spallation products 
that reconstruct inside the analysis cylinder to the number that 
reconstruct in the full LS volume. This fraction combined 
with the total LS volume in the balloon, which was measured 
directly during construction to be 1171 ±25 m 3 , gives a fidu- 
cial volume of 176.4m 3 with an uncertainty of 3.1%. When 
all selection cuts, to be described, are applied the remaining 
exposure is 123 kton-days. 

This paper focuses on neutrino-electron elastic scattering, 
v + er — > v + e~ (ES) due to 8 B solar neutrinos. With 
3.423 x 10 32 e~ targets per kton of LS, we expect, including 
oscillation, 1.19 recoil events per kton-day from 8 B in the 
energy-region of interest (ROI) 5.5-20 MeV. In our calcula- 
tion we adopt the standard solar model (SSM) of Serenelli et 
al. Il6ll which uses the solar abundances of Asplund et al. IU7I1 
(AGSS09). This predicts a total 8 B neutrino flux, independent 
of neutrino flavor, of $s B 



4.85±g;|| xlO 6 * 



the ES cross section from [18]; the B-i> e spectrum of Winter 
et al. lfl9tl : the oscillation parameters from the global analy- 
sis analysis in 01; and fold in the detector response. ES from 
hep neutrinos and neutrino interactions on carbon in the LS 
are treated as a background and are discussed later in the pa- 
per. We note that the best choice of solar abundances to use 
in the SSM is an unresolved question. In AGSS09 the he- 
lioseismic discrepancy which emerged when the abundances 
of Asplund et al. ll2oTl (AGS05) were adopted still persists, al- 
though it is not as severe. If the older, 1 -dimensional, solar 



abundance evaluation of Grevesse and Sauval ]21|](GS98) is 
used, the SSM has excellent consistency with helioseismology 



and the predicted 8 B flux is <&s B = 5. 



-0.65 



xl0 b cm- 



Assuming no sterile neutrinos, the total 8 B flux is directly 
observed in the SNO neutral current measurement to be 
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KamLAND data 
Best Fit 
Bv ES 
J2 "Be decay 
Li decay 
External y's 
'C(y e , e) 13 N 
□ Other 



TABLE TIT: The systematic uncertainties associated with the un- 
binned fit to the energy spectrum of the 8 B candidates. The detection 
efficiency is dominated by our fiducial volume uncertainty. 



11 12 13 20 
Reconstructed E [MeV] 



FIG. 1 : Energy spectrum of 8 B candidates with the best-fit spectrum 
and background components from the unbinned energy and rate anal- 
ysis. The histograms display the results in bins of 0.5 MeV except for 
the last bin which due to limited statistics extends from 13.5 MeV to 
20 MeV. 



$8 B = 5.140i^207xl0 6 cm- 2 s- 1 fToh . 

With the 5.5 MeV analysis threshold the background is 
dominated by decays of light isotopes produced by muon spal- 
lation. An in-depth study of muon activation at KamLAND 
can be found in 1U5I1 . A study of light isotope production 
shows that most (> 80%) light isotope backgrounds are cor- 
related with muons which produce more than 700,000 photo- 
electrons (p.e.) in the PMT array. We denote these as bright 
muons. The rate of bright muons is 0.037 Hz, while the total 
rate of muons passing through the LS is 0.198±0.014Hz. We 
apply a number of muon-related cuts to reduce these spallation 
backgrounds. All events within 200 ms of a preceding muon 
are rejected. This time veto of the full detector significantly 
reduces the background from 12 B/ 12 N. For non-bright muons, 
for which the muon tracking algorithm converged successfully 
— well-tracked muons — a 5 s veto is applied within a 3-m- 
radius cylinder around the muon track. Using the 12 B/ 12 N 
candidates we determine that 97±6% of spallation products 
are contained within this cylinder. For bright muons and any 
muon with a poorly reconstructed track, the 5 s veto is ap- 
plied to the full detector. These cuts reduce the exposure by 
62.4±0.1%, to a total exposure of 123 kton-days. The spal- 
lation background events remaining after these cuts are ex- 
pected to come mainly from long lived (r > Is) spallation 
products, 11 Be, 8 Li, 8 B. The total production rates of these 



key isotopes 1 15], along with the rates correlated to non-bright 
muons, are summarized in Table U 

The next-largest background is from external gamma rays 
which are primarily the result of (71,7) reactions in the sur- 
rounding rock cavern and stainless steel detector elements. 
The external 7-ray spectrum has peaks in reconstructed en- 
ergy at 8.5 MeV and 10 MeV from stainless steel, and at 
5.5 MeV from neutron capture on silicon in the rock. The 
cylindrical fiducial volume was chosen to optimize the shield- 
ing for a given exposure. The closest points to the cylindri- 



Source 



Uncertainty (%) 



n Be 
8 Li and 8 B 
External gamma rays 
Other Backgrounds 
Detection Efficiency 
Energy Scale 



10.8 

3.3 

6.8 

1.1 

6.3 

0.8 



cal external rock cavity are at the balloon equator, while the 
closest stainless steel component is a balloon support struc- 
ture at the top of the detector. A GEANT4-based Monte 
Carlo Ir22ll23ll . including the full detector and shielding geom- 
etry, simulated the effect of LS self-shielding from sources in 
the stainless steel and the rock. The simulation indicates that 
gammas are attenuated approximately exponentially with an 
attenuation length of 53.0±0.1 cm for rock and 50.7±0.1 cm 
for stainless steel. Using data within a 6-m-radius volume, 
we observe attenuation lengths for gammas from the rock and 
stainless steel of 59. Oil. 9 cm and 54. 9± 1.9 cm respectively. 
We estimate 25. 2± 12.6 electron-recoil-like events from exter- 
nal 7-rays in the ROI within the cylindrical fiducial volume. 
The uncertainty in the estimate comes from the difference in 
the observed and simulated attenuation lengths. 

As remarked earlier, events from hep solar neutrinos and 
solar neutrino interactions on carbon are treated as a back- 
ground. Using the SSM with AGSS09 and the hep spectrum 
from 124J, we estimate, including oscillation, 0.6± 0.1 elec- 
tron recoil events from hep neutrinos in the ROI. In our calcu- 
lation we use the ES cross section, neutrino oscillation param- 
eters and detector response as was done for the 8 B ES calcu- 
lation. The uncertainty is dominated by the difference in the 
flux prediction of the SSM with the AGSS09 versus with the 
GS98 solar abundances. 

There are 4.30xl0 31 carbon nuclei per kton of LS if 
we assume a natural 13 C of 1.10%. Using the cross sec- 
tions calculated in I25I1 . we find the largest is e -C scattering 
background contribution to be from charged current (CC) 
scattering, 13 C+^ e — > 13 N+e~ from 8 B. We estimate, includ- 
ing oscillation, that scattering to the ground state of 13 N 
produces 5. 8 ±1.4 events in the ROI; and scattering to the 
3.51 MeV first excited, which decays by proton emission, con- 
tributes 1 . 1 ± 0.4 electron-recoil-like events in the ROI. In this 
estimate an additional uncertainty of 30% on the cross sec- 
tion is included. The contribution from higher states of 13 N, 
neutral current (NC) scattering by 8 B-^ and hep — v NC and 
CC interactions on carbon is estimated to be less than 0.13 
events in the ROI, and is considered negligible given the other 
backgrounds. These calculations use the same solar-z/ fluxes, 
spectra, detector response and oscillation parameters as be- 
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FIG. 2: Summary of measurements of the 8 B flux using neutrino-electron scattering: this work KamLAND, Kamiokande II flUIT^I. Super 
Kamiokande I (J, Super Kamiokande II @], Super Kamiokande III (H, SNO D 2 HE3], SNO Salt S[0i, SNO NCD GJ, and Borexino 



fl4ll . The SNO neutral current measurements ll9l lToll are shown for reference with closed triangles. 



fore. The effect of using the oscillation parameters from the 
global analyses in 113112611 were investigated but do not change 
the result due to the large flux and cross section uncertainties. 

In addition to these neutrino interactions, we expect a small 
background from atmospheric neutrino interactions. Assum- 
ing a flat spectrum, we use the candidate events in the range 
of 20 MeV to 35 MeV to extrapolate the atmospheric neutrino 
contribution in the signal region, and find 2.0±2.0 events. 

A background from positrons from inverse beta decay in- 
duced by reactor antineutrinos is tagged by coincidence with 
delayed neutron capture 7-rays. The mean capture time in 
KamLAND is -~207 /is. The tagging cuts are: a delayed en- 
ergy cut, 2.04MeV<E<2.82MeV, to select the 2.2MeV n-p 
capture 7; a timing cut, 0.5/is<AT<660/is; and a position 
cut, AR< 1 .6m. AT and AR refer respectively to the time and 
distance between the prompt positron and delayed capture 7. 
The tagging efficiency is 0.940±0.006. The reactor antineu- 
trino flux at KamLAND is dominated by Japanese power re- 
actors, we calculate the flux and spectrum using data provided 
by the Japanese power companies and a simplified reactor 
model |27l| . The residual background is 1.6 ±0.1, where the 
uncertainty in the neutrino oscillation parameters is included. 

The total background, broken down in Table ITT1 is estimated 
to be 157.3±23.6 events. While 8 Li has the largest produc- 
tion rate, 11 Be, due to its long half-life (ti 12 = 13.8 s), is the 
largest contribution to the total background. 

After applying all the cuts, 299 electron recoil candi- 
dates are found in the ROI in the 123 kton-day exposure. 
Subtracting the 157.3±23.6 background events gives a 8 B 
rate of 1.16±0.14(stat)±0.21(syst)events per kton-day which 
is consistent with our prediction including neutrino oscil- 
lation. If we neglect neutrino oscillation, assume a pure 
v e flux, and correct for the 5.5 MeV threshold, the mea- 
sured rate corresponds to a spectrum integrated flux of 
$ ES = 2.17±0.26(stat)±0.39(syst) x 10 6 cm^s" 1 . 

A second analysis of this data uses an unbinned maximum 
likelihood fit to the 8 B candidate energy spectrum. The nor- 



malization of the individual backgrounds are allowed to vary 
but are constrained by a penalty to the expected value. The 
best-fit rate is 1.49±0.14(stat)±0.17(syst)events per kton- 
day, with a goodness-of-fit of 49%. Once again assuming a 
pure v e flux, this corresponds to a spectrum integrated flux 
of $ ES =2.77±0.26(stat)±0.32(syst) x 10 6 cm^s" 1 . Fig. [Q 
shows the energy spectrum of the candidate events with the 
best-fit solar neutrino and background spectra. This anal- 
ysis does not include the small change in the shape of 
the neutrino spectra due to oscillation. Including this ef- 
fect does not significantly change the best-fit rate, <&es = 
2.74±0.26(stat)±0.32(syst)xl0 6 cm~ 2 s- 1 , with a goodness 
of fit of 57%. 

The KamLAND result is compared to other measurements 
of the total 8 B flux through neutrino-electron elastic scat- 
tering in Fig. [2] For this comparison the reported statis- 
tical and systematic errors are added in quadrature. The 
mean of the experiments, weighted by their uncertainties, is 
<I> E s =2.33±0.05 x 10 6 cm" 2 s" 1 . This is dominated by the 
very precise Super-Kamiokande water-Cherenkov measure- 
ment. The reduced x 2 /NDF is 6.6/ 8 where the KamLAND 
Rate + Energy result and the Borexino E>3MeV result are 
used in the calculation. 

A key prediction of the LMA-MSW solution to the so- 
lar neutrino problem is the expected transition from matter- 
dominated to vacuum-dominated oscillations, depending on 
the neutrino energy and the region of the solar interior probed 
by the neutrinos as they journey out of the sun. A general fea- 
ture is that the survival probability, defined as the ratio of the 
oscillated j/ e -flux to the unoscillated flux, tends to increase 
with decreasing solar neutrino energy. This effect is the re- 
sult of neutrino species ( 8 B, 7 Be, pp, etc.) either moving off, 
or being produced in a region not satisfying the MSW reso- 
nance condition. These neutrinos are then described by vac- 
uum oscillation. The curve in Fig. [3] shows the ratio of the 
expected neutrino-electron elastic scattering rates with and 
without oscillation calculated using the AGSS09-SSM. The 
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and Japan Atomic Energy Agency. The Kamioka Mining and 
Smelting Company has provided service for activities in the 
mine. 
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Electron Recoil Energy [MeV] 



FIG. 3: Ratio of the measured (threshold corrected) solar neutrino- 
electron elastic scattering rate and that predicted by the AGSS09- 
SSM as a function of energy (data points). The average of the water- 
Cherenkov detector results above 5 MeV is shown. The horizontal 
error bars indicate the range of energy over which the flux was mea- 
sured. The central value indicates the flux averaged energy for that 
range. The solid curve (and ± 1 a band) shows the ratio of AGSS09 
with neutrino oscillation included to the prediction without oscilla- 
tion. The discontinuities in the prediction arise as different neutrino 
sources become significant, the large uncertainty in the CNO neutri- 
nos is evident at ~ 1.2 MeV. 



data points are the survival probabilities deduced, relative to 
the SSM prediction, from neutrino-electron elastic scatter- 
ing measurements at Borexino, KamLAND and the combined 
water-Cherenkov experiments. The Borexino 8 B measure- 
ment with a 3.0 MeV threshold Il4ll and the 7 Be flux mea- 



surement 1281 12911 are consistent with LMA-MS W prediction. 
The measurements are not yet precise enough to resolve the 
issue of solar abundances. 

This letter reports a measurement of the 8 B solar neu- 
trino flux measured through neutrino-electron elastic scatter- 
ing in the KamLAND LS. Our result is in good agreement 
with existing measurements from both water-Cherenkov and 
liquid-scintillator-based detectors. In the data-set presented, 
208 T1 background limits the KamLAND analysis threshold to 
5.5 MeV. We expect to achieve a lower threshold in a future 
analysis of data taken after purification of the KamLAND LS. 
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